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Abstract - In France, implementation of systematic screening programs in 2000. as a complement 
to the mandatory reporting of animals with clinical signs of BSE (passive surveillance), revealed 
certain limitations of the mandatory system. Indeed, systematic screening showed that some BSE 
cases were not detected by the clinical surveillance system, implying considerable BSE case under- 
reporting throughout the epidemic. As the most likely explanation for variant Creutzfeldt-Jakob 
disease (vCJD) is exposure to the aetiologic agent of BSE, it is essential to reconstruct the French 
BSE epidemic pattern accounting for this under- reporting. We estimated age- and year-specific 
incidence rates of BSE by using a back-calculation method. This approach relies on the principle 
that the number of clinical BSE cases is the consequence of the number of BSE-infected animals 
after a known incubation time, defined as the time between infection and clinical onset. We 
generalized this model to take into account epidemiological characteristics of BSE, such as French 
cattle mortality, BSE case reporting probability, and age-dependent susceptibility and/or exposure 
to the BSE agent. We confirmed that the average BSE incubation period is five years and that the 
peak risk of bovine infection occurs between 6 and 12 months of age. The results also showed that 
the proportion of underreporting is the most influential parameter in the model, and that BSE was 
substantially underreported until rapid tests were introduced. Indeed, only 103 BSE cases were 
detected by passive surveillance up to June 2000, while we estimated that there was 301 200 (95% 
confidence interval (CI) [27 600-837 600]) cattle infected by the BSE agent. Despite uncertainty 
over the beginning of the epidemic, we showed that the French BSE epidemic in the late 1980s was 
completely undetected, and only the second wave, after 1990, was observed. 

epidemiology / BSE / under-reporting / vCJD / back-calculation 



1. INTRODUCTION 2001, the surveillance system underwent 

several changes. From mid-2000, in addi- 

In France, BSE became a notifiable dis- tion to the mandatory reporting system, a 

ease in June 1990. In December 1990, a pilot study of rapid testing was implemented 

mandatory passive surveillance system was on cattle at risk (dead-on-farm cattle, emer- 

set up, in which veterinary practitioners and gency slaughtered cattle and euthanatized 

farmers were required to report animals cattle). Then, from January 2001, system- 

with clinical signs. The first case of BSE atic screening was extended to all cattle 

identified by this surveillance system was over 30 months of age entering the food 

detected in 1991. Between 1991 and June chain, and this age was reduced to 24 months 

2000, 103 cases were identified by passive in July 2001. Implementation of active sur- 

surveillance [2]. From mid-2000 until July veillance revealed the inefficiency of the 
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passive surveillance, both in France and 
elsewhere. In France, comparison of the pas- 
sive and active surveillance systems between 
July 2001 and June 2002 showed that only 
20% of BSE cases were identified by pas- 
sive surveillance. This implied considerable 
BSE case under-reporting throughout the 
epidemic. 

It is now widely recognised that the most 
likely explanation for human variant Creut- 
zfeldt-Jakob disease (vCJD) is exposure to 
the BSE agent [10, 26, 33]. As of 5th 
January 2004, there were 145 cases of def- 
inite or probable vCJD [16] in the United 
Kingdom. The future size of the vCJD epi- 
demic in the United Kingdom is specula- 
tive, although the most recent predictions 
[25, 27, 31] based on statistical modelling 
are more optimistic than previously. The 
number of vCJD cases in France - 6 con- 
firmed cases [28] - is far too small to 
attempt direct modelling of the vCJD epi- 
demic in the French population. Neverthe- 
less, a comparative assessment of the tem- 
poral pattern of exposure to the BSE agent 
in the United Kingdom and France could 
help to estimate the size of the vCJD epi- 
demic in France. The French population has 
been exposed to the BSE agent from three 
sources: via travel to the United Kingdom 
and British beef imports since 1980, the 
likely start date of the BSE epidemic [13], 
and via the BSE epidemic in French cattle. 
A study of French blood donors showed that 
travel to the United Kingdom contributed 
very little to the overall risk in the French 
population [5]. Exposure of the French 
human population via imported British 
bovine products was estimated to represent 
between 5% and 10% of the exposure level 
of the British population [4]. This exposure 
source only posed a high risk during the 
period 1980-1989, i.e. the period before the 
British ban on specified offal from all cattle 
entering the human food chain (November 
1989). In France, measures intended to pre- 
vent bovine infection by the BSE agent 
were taken from 1989. However, it was only 
in June 1996 that high-risk bovine tissues 
("specified offal") were excluded from 



human consumption. Thus, the French pop- 
ulation was exposed to the BSE agent via 
French bovine products until June 1996. 
However, only animals slaughtered for con- 
sumption late in the incubation period are 
likely to be infectious. 

Another key parameter is necessary to 
study the exposure of the French popula- 
tion, namely the start date of the French 
BSE epidemic. BSE was first recognised in 
England in 1985. Supplementary feeds con- 
taining meat and bone-meal (MBM) were 
identified as a potential vector of the BSE 
agent. A French Senate report [8] states that 
MBM was already being imported by France 
from the United Kingdom in 1985 (no infor- 
mation is available before this date) and, 
according to a European report, from other 
European countries as early as 1980 [30]. 
As the earliest birth cohort of a BSE case 
dates back to 1983, French cattle were 
already potentially exposed to the BSE 
agent at that time. For the purposes of this 
study, we assumed that the French epidemic 
began in 1980; very similar results were 
obtained when a start date of 1 983 was used. 

We used the back-calculation method 
[9] to estimate longitudinal trends in the 
incidence of BSE infection in France, and 
to deduce the number of infectious animals 
slaughtered for human consumption. This 
method relies on the principle that the 
number of clinical BSE cases is the conse- 
quence of the number of BSE-infected ani- 
mals after a known incubation time, defined 
as the time between infection and clinical 
onset. We generalised this model to take 
into account epidemiological characteris- 
tics of BSE, such as a time-dependent BSE 
case-reporting function. 



2. MATERIALS AND METHODS 



2.1. Data 

Our analysis was restricted to the 103 bovine 
clinical cases [2] detected by passive sur- 
veillance up to June 2000 (Tab. I). All the 
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Table I. Distribution of clinical cases, detected by passive surveillance between January 1991 and 
June 2000, by age class and time at clinical onset. 



^"^-^^^ Age class 

^^^^^^^(years) 

Period ^\^^^^ 
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Jan.-June 1997 



















1 











1 








2 


July-Dec. 1997 
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Jan.-June 1998 
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1 
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July-Dec. 1998 
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1 
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13 


Jan.-June 1999 
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3 


4 


3 























12 


July-Dec. 1999 
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7 


2 


1 
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1 





1 





19 


Jan.-June 2000 
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4 


6 


4 


2 


1 





1 








23 


Total cases by age 


: class 


1 


13 


17 


19 


17 


9 


6 


5 


5 


4 


4 


3 


103 



Jan.: January; Dec: December. 



animals were cows aged between 3.9 and 
9.4 years at clinical onset; 85% were aged 
between 4 and 7 years, and mean age at clin- 
ical onset was 5.9 years. 

Specific developments would be needed 
to take into account clinical cases detected 
after June 2000, as the multiple changes to 
the surveillance system between July 2000 
and July 200 1 could have introduced detec- 
tion biases. Indeed, a retrospective clinical 
investigation [12] showed that some posi- 
tive animals found in the pilot study should 
have been included in the mandatory report- 
ing system because they had clinical signs 
of BSE. In addition, active surveillance 
improved awareness of BSE and the effi- 
ciency of mandatory reporting [14]. A pos- 



sible "escape route" for BSE cases was the 
over 30-month-cattle destruction program. 
Indeed, between January 200 1 and June 200 1 , 
more than 179 000 cattle over 30 months of 
age were destroyed without being tested. 
No simple models can account for such 
rapid changes, and we therefore decided 
first to develop a model based only on the 
period up to June 2000. 

2.2. BSE epidemiological 
characteristics 

Epidemiological data on clinical onset 
by age (Tab. I), assuming a long BSE incu- 
bation period [6, 24, 32], suggest that most 
infections occur shortly after birth. We 
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mean=5.0 years 
variance=1 .Byears^ 




Figure 1. BSE epidemiological input data. (A) Estimated survival distributions of French (:A-) and 
British (-Q-) cattle (23). (B) Estimated probability density function of the BSE incubation period 
(q = 1.667; p = 5). (C) Estimated cumulative distribution function of cattle age at infection, in six 
classes (in years) ]0-0.5], ]0.5-l], ]l-2], ]2-3], ]3-5], and ]5-30]. (D) Time-dependent BSE case 
reporting function, assuming a 5% reporting probability in June 2000. In the sensitivity analysis, we 
considered p-values from (-•-) to -1 ( — ); p = -0.6 (-■-) minimised the model-selection criterion. 



therefore used age-dependent susceptibil- 
ity and/or exposure to infection, as availa- 
ble data did not allow us to discriminate 
between the two. 

In France, 49% of cattle are slaughtered 
before age one year (Fig . 1 A, ^) . Thus, age 
at slaughter and natural mortality both need 
to be incorporated, as this might censor 
clinical onset in infected animals, or some 
animals could be slaughtered before being 
exposed to the infectious agent. 

The sharp increase in BSE incidence 
when the first systematic screening pro- 
gram was instigated revealed the poor effi- 
ciency of passive surveillance; a time- 
dependent BSE case reporting function 
must thus be used. As only 20% of all BSE 
cases were identified by passive surveil- 
lance between July 200 1 and June 2002, we 



used an upper limit of 20% for the reporting 
probability in June 2000. 

2.3. Generalised back-calculation 
method 

The back-calculation method was extended 
to take into account the effect of age, using the 
method described by Becker and Marshner 
[7, 17], the cattle survival distribution, and 
a time-dependentprobability of BSE report- 
ing. 

Let Na^f and Y^ , be the random numbers 
of new infected animals and new cases with 
age a at time t, and E{N^,), EiY^f), their 
expectations. Assuming that the distribu- 
tion incubation time has a density /(?), that 
S(a/a') represents the probability that an 
animal survives to age a from all causes of 



Case 1:06-cv-00544-JR Document 14-1 1 Filed 11/03/2006 Page6of15 



The unrecognised French BSE epidemic 



353 



mortality, knowing that it was alive at age 
a ' (age at infection) and that A (f) is a time- 
dependent probability that a given BSE 
case is actually reported at time t, the model 
becomes: 






,)f{t-s)S{a/a-t + s)dsA{t). 



The distribution incubation time was 
supposed to be independent on infection 
age and the survival distribution was sup- 
posed to be independent on time. 

The incubation period was assumed to be 
gamma distributed, so the probability density 
function of the incubation period/(r) was: 



m = 



Vf<2 
l /(f-2f"^xp(-g(^-2)) 

I ro,) 



Vf>2. 



In order to reproduce the observed time- 
lag of two years before any clinical cases 
are seen, we assumed that the probability of 
a B SE incubation period of two years or less 
was nil. The distribution mean was {p/q+2) 
years and the distribution variance was 
p/q2. 

The survival distribution (Fig. 1 A, ^) was 
estimated from three data sources: Office 
National Interprofessionnel des Viandes, 
de I'Elevage et de 1' Aviculture (OFTVAL), 
Agence Fran9aise de Securite Sanitaire de 
I'Alimentation (AFSSA) and Direction 
Generale de I'Alimentation (DGAL). The 
OFTVAL report [29] gave the number of 
calves (aged between and 6 months), and 
the number of young cattle (aged between 
6 and 24 months), slaughtered in abattoirs in 
2001. The DGAL dataset [18] gave the 
number of calves and young cattle sent for 
rendering in 2000. These last two data 
sources allowed us to estimate survival up 
to 24 months. The AFSSA dataset [1] gave 
the number of cattle aged over 24 months 
which were sent to abattoirs or rendering in 



2001, allowing us to estimate the survival 
of cattle aged over 24 months. We derived 
conditional survival probabilities from this 
survival distribution. We considered sur- 
vival given being alive at the time of infec- 
tion, as most cattle die or are slaughtered 
before age one year in France. Thus, the prob- 
ability of being alive at five years of age given 
being alive at one year age is not equal to 
the probability of being alive at five years age. 
A logistic form was assumed for the 
time-dependent reporting function, starting 
in June 1990, with zero reporting before this 
date as there was no monitoring system: 



A(0 = 

[ exp(e + [i(Ji/«e2000-f)) 
\l + exp{e + ^{June 2000 - 1)) 





Vf>1990, 
Vf<1990. 



The logistic curve depended on two 
parameters, p the form parameter and 9 the 
parameter determining the reporting prob- 
ability in June 2000. Let A (/MMe2000) be 
the probability of being reported in June 
2000, then 



= In 



A(/Mwe2000) 
vl -A(June2000 



The unknown time- and age-specific 
BSE infection numbers were modelled 
using the multiplicative model E{N^ g) = 
aj^s' these parameters are defined except 
for a multiplicative constant; X^ was the 
time-varying risk of feed-borne infection 
and a^ was the age-dependent susceptibil- 
ity/exposure. However, we did not con- 
sider age-specific susceptibility/exposure but 
rather susceptibility/exposure per age class. 
In the absence of epidemiological data to 
suggest a parametric family, we did not 
specify a parametric form for the time-var- 
ying BSE-infection risk Xy 

Under the assumption that the age- and 
time-specific new BSE-infected animal num- 
bers, N^ p were independent Poisson vari- 
ates then the age- and time-specific new 
clinical case numbers 7^ , were independent 
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Poisson variates. Then, we had the likeli- 
hood function corresponding to the age- 
and time-specific observed clinical BSE 
cases yai- 

L=Lia,X\y) = hn'^^^^^^^=^ 

■'^ at 

where n^_, = E{Y^;). 

Non parametric maximum likelihood 
estimation to A,j was implemented using the 
EM algorithm [15] to which we added a 
smoothing step (EMS algorithm [7]). The 
time-varying BSE-infection risk, X^, was esti- 
mated annually from 1980 to 1996. Years 
were defined so that, for example, 1980 
consists of the period between 1 July 1980 
and 30 June 1981. As we considered clini- 
cal BSE cases detected up to June 2000, we 
could not estimate the time-varying BSE- 
infection risk, or BSE infection numbers, 
beyond June 1997. Indeed, because of the 
long BSE incubation period there is little 
information contained in BSE incidence 
data about the numbers of animals infected 
most recently. The continuous estimates 
were obtained from the discrete estimates 
by linear interpolation. 

The back-calculation method allows one 
to estimate the past time- and age-specific 
BSE infection number, ZifA'^^j, from clini- 
cal BSE case data, £(^0;), provided that 
other model parameters are known. How- 
ever, the susceptibility/exposure at infec- 
tion was unknown and the time-dependent 
reporting function and the incubation period 
density depended on unknown parameters. 
So we performed systematic sensitivity 
analyses for each unknown parameter. 

2.4. Sensitivity analyses 

We explored the sensitivity of the model 
results according to unknown parameter 
values. 

We varied the mean incubation period 
distribution between 4 and 6 years and the 
variance between 1 and 3 years^. 



We fitted a flexible function to the time- 
dependent reporting function, dependent 
on two parameters. We considered values 
between and -1 for the form parameter, 
p. This made it possible to explore a wide 
range of reporting forms, ranging from con- 
stant reporting throughout the epidemic 
(P = 0) to a reporting function with very low 
reporting probabilities in the early stages of 
the epidemic and a sudden increase in the 
recent past (p = -1). By varying the second 
parameter, 9, we varied the reporting prob- 
abihty in June 2000 between 5% and 20%. 

A representative range of susceptibility/ 
exposure forms was examined, such as con- 
stant exposure susceptibility/exposure, sus- 
ceptibility/exposure decreasing with age, 
and no susceptibility/exposure except in the 
first six month of life. Also, we varied the 
number and spacing of the age class knots. 
We considered six-class susceptibility/ 
exposure (in years): ]0-0.5], JO.5-1], ] 1-2], 
]2-3], ]3-5], J5-30], and compared them to 
five-class susceptibility/exposure: ]0-l], 
] 1-2], ]2-3], ]3-5], ]5-30], and seven-class 
susceptibility/exposure: JO-0.5], ]0.5-l], 
]l-2], ]2-3], ]3-5], ]5-7], ]7-30]. 

The best model was selected by using 
Akaike's Information Criterion (AIC) [11]: 
AlC(model) = -2 x log(L) H- 2 x K where 
L is the model's likelihood and K is the 
number of estimated parameters. We con- 
sidered the empiric rule of Burnham and 
Anderson [11], which retains a model if 
(AlC(model) - MinAIC) < 2 where MinAIC 
is the minimum AIC, i.e. the AIC of the best 
model. 

All unknown parameters were selected 
to minimize the AIC. However, as all the 
reporting probabilities in June 2000 that we 
explored were compatible with a good fit 
(similar AIC), we needed external data to 
select among the hypotheses. Therefore, to 
establish the reporting probability in June 
2000, so to determine the parameter 9, we 
compared annual predicted clinical BSE 
cases with observed clinical cases. 

To determine precision, we used boot- 
strap [21] estimates considering that the 
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observed age- and year-specific incidence 
of BSE cases arose from a Poisson distri- 
bution. The 95% confidence interval for 
each parameter was obtained by excluding 
the 5% most extreme values obtained from 
100 bootstrap samples. 

2.5. Predictions 

To obtain the number of new clinical 
cases at time t, we used the following func- 
tion: 

A 

\E(YJ(h = 



!! 



E(N^_,^^^)f(t-s)S{a/a-t + s)dsK(t)da. 



Up to June 1997 we used estimates of 
EiN^ ,) and beyond June 1997 we assumed 
that no new infections occurred. Regarding 
the time-dependent reporting function, beyond 
June 2000 we assumed that the reporting 
probability of clinical cases was 20%. 

2.6. Number of late-stage BSE cases 
slauglitered for consumption 

To compare exposure of the British and 
French populations to the BSE agent, we 
defined late-stage animals as animals slaugh- 
tered within 1 2 months of clinical onset, as 
in the British study [23]. From the estimated 
number of BSE-infected cattle, we ascer- 
tained by simulation the number of late- 
stage animals slaughtered for consumption. 
To each infected animal, we randomly 
assigned an age at infection, an incubation 
period, and a lifetime, given that the animal 
was alive at the age of infection. Bearing in 
mind mortality directly attributed to clini- 
cal cases, we imposed that lifetime ranges 
between age at infection and age at infec- 
tion plus incubation period. 

We thus obtained the number of BSE- 
infected cattle between times s and s + A 
aged between a and a H- A at infection with 



an incubation period between d and d + A, 
and a lifespan between a' and a'+ A. We 
were thus able to deduce the number of 
infected cattle slaughtered between times t 
and t + A{t = s + {a' - a)), between x and 
X + A times before the clinical onset of BSE 
(x = d — (a' - a)), and aged between a ' and 
a'H- A at slaughter. 

The method of back-calculation, the EMS 
algorithm and the simulations were imple- 
mented using computer programs written in 
C++ for this purpose. 

As the likelihood is no longer maximised 
when the smoothing step is included in the 
EM algorithm, we used a convergence cri- 
terion based on the values of the parameter 
as Becker and Marshner [7]. 



3. RESULTS 

The best model, according to Akaike's 
Information Criterion, suggests that the 
average BSE incubation period is five years 
(Fig. IB) and that most infections occur 
between 6 and 12 months of age (Fig. IC). 

With reference to the form parameter of 
the reporting curve, p = -0.6 minimises the 
AIC, we also present results for (3 = -0 .5 and 
P = -0.7 as they satisfy the empiric rule of 
Burnham and Anderson. To establish the 
reporting probability in June 2000, we com- 
pared annual predicted clinical BSE cases 
(Tabs. II and III) with observed clinical cases. 
Passive surveillance detected 134 clinical 
BSE cases between July 2000 and June 
2001 and 60 clinical BSE cases between 
July 2001 and June 2002 (respectively 75 
and 161 BSE cases were detected among 
cattle at risk). Table II shows short-term pro- 
jections of the clinical BSE case incidence 
according to the form parameter and the 
level of the reporting curve in June 2000. 
Table III illustrates longer-term projections 
according to the reporting probability in 
June 2000, the form parameter p being fixed 
at -0.6. Comparison of predicted cases with 
the 134 and 60 clinical cases detected in the 
periods from July 2000 to June 2001 and 
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Table II. Predicted incidence of clinical BSE cases in French cattle (with 95% bootstrap CIs) for the 
periods July 2000 to June 2001 (2000-2001) and July 2001 to June 2002 (2001-2002) according to 
the form parameter, P, and the reporting probability in June 2000 (A(June2000) ). We assumed that 
no animals born after June 1997 became infected. 



Predicted incidence of clinical BSE cases (with < 



,CI) 



|3 = -0.7 



P = -0.6 



P = -0.5 







2000-2001 


2001-2002 


2000-2001 


2001-2002 


2000-2001 


2001-2002 


A(Jiine2000) 


= 0.05 


139 


69 


132 


67 


127 


64 






[99-207] 


[47-108] 


[95-199] 


[45-105] 


[90-191] 


[43-102] 


A(Jiine2000) 


= 0.10 


69 


34 


66 


33 


63 


32 






[49-103] 


[23-54] 


[47-99] 


[22-52] 


[45-95] 


[21-51] 


A(/w«e2000) 


= 0.15 


45 


23 


43 


22 


41 


21 






[32-68] 


[15-35] 


[31-65] 


[15-34] 


[29-63] 


[14-34] 


A(/w«e2000) 


= 0.20 


33 


17 


32 


16 


31 


16 






[24-50] 


[11-26] 


[23^8] 


[11-26] 


[22-47] 


[11-25] 



Table III. Predicted annual incidence of asymptomatic and clinical BSE cases, (with 95% bootstrap 
CIs) according to the reporting probability in June 2000 (A{June2000y). We assumed that parameter 
p of the reporting curve was -0.6 and that no animals born after June 1997 became infected. 

Predicted annual incidence of asymptomatic and clinical BSE cases (with 95% CI) 
A(/u«e2000) =0.05 A(June2000) =0.10 A(Jiine2000) =0.15 A(JMwe2000) =0.20 



July 2002/ 


127 [82-205] 


63 [40-102] 


42 [27-68] 


31 [20-50] 


June 2003 










July 2003/ 


40 [25-66] 


20 [12-33] 


13 [8-22] 


10[6-16] 


June 2004 










July 2004/ 


11 [6-18] 


5 [3-9] 


4 [2-6] 


3 [2-4.] 


June 2005 










July 2005/ 


2 [1-4] 


1 [1-2] 


1 [0-1] 


1[0-1] 


June 2006 










July 2006/ 


1 [0-1] 


[0-0] 


[0-0] 


0[0-0] 


June 2007 











July 200 1 to June 2002, respectively, reveals 
that the reporting probability in June 2000 
was no higher than 10% (Tab. II). In addi- 
tion, the comparison of annual predicted clin- 
ical or asymptomatic cases for the period 
July 2002 to June 2003 (Tab. Ill) with the 
172 clinical or asymptomatic BSE cases 
recorded between July 2002 and June 2003 
shows that the assumption of a reporting 



probability of 10% or more in June 2000 
would underestimate the true situation. 
Thereafter, we inferred that the best esti- 
mate of the reporting probability in June 
2000 was 5% and that the form parameter of 
the reporting curve was -0.6 (Fig. ID, -■-). 

Under the previous assumptions, we 
estimated that 301200 (CI = [27 600- 
837 600]) cattle were infected by the BSE 
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Table IV. Estimated numbers of infections in French cattle (with 95% bootstrap CI) from June 1980 
to June 1997 (A) and from June 1987 to June 1997 (B) according to the form parameter, P, and the 
reporting probability in June 2000 (A{June2000)). 

Estimated numbers of infections in French cattle (with 95% CI) 



P = -0.7 



P = -0.6 



P = -0.5 



(A) 

A{Jiine2000) =0.05 

A{Jiine2000) =0.10 

A{Jime2000) =0.15 

A{Jime2000) =0.20 

(B) 

A{Jiine2000) =0.05 

A(Jime2000) =0.10 

A{Jime2000) =0.15 

A(JMne2000) =0.20 



656 400 
[41600-1872 800] 

311400 
[20 000-888 000] 

196 500 
[12 700-560 000] 

138 900 
[9100-395 700] 

78 900 
[34 500-135 400] 

37 700 
[16 600-64 700] 

24 000 
[10 600-41 100] 

17 100 
[7 600-29 300] 



301 200 
[27 600-837 600] 

143 000 
[13 300-397 300] 

90 300 
[8 500-250 800] 

64 000 
[6100-177 300] 

51300 
[24 300-84 700] 

24 600 
[11700^0 600] 

15 700 
[7 500-25 900] 

11300 
[5 400-18 500] 



141 100 
[18 800-377 200] 

67 200 
[9100-179 200] 

42 500 
[5 800-113 200] 

30 200 
[4 200-80 200] 

34 200 
[17 300-54 800] 

16 500 
[8 400-26 400] 

10 600 
[5 400-16 900] 

7 600 
[3 900-12 200] 



agent since 1980 (Tab. IV(A)). Table IV 
shows estimates of BSE-infected cattle 
numbers according to the two parameters of 
the reporting function. It should be noted 
that these estimates are strongly dependent 
on reporting probabilities. Also, the under- 
reporting proportion was the parameter of 
the model with most influence on these esti- 
mates. 

The confidence intervals (Fig. 2 A, dot- 
ted lines) reflect the uncertainty in the esti- 
mated number of infections occurring dur- 
ing the 1980s. From 1987 onwards, estimates 
become more precise as the confidence 
interval amplitude is lower. The infection 
number rises between 1987 and 1990, then 
falls between 1990 and 1992; 1992 to 1995 
saw another rise but smaller than the previ- 
ous one and, finally, we observed a new fall 
after 1995. Between July 1987 and June 



1997, an estimated 51 300 (CI = [24 300- 
84 700]) cattle were infected (Tab. IV(B)). 
Very similar estimates were obtained when 
the epidemic start date was assumed to be 
1987 (results not shown). 

From these infection number estimates, 
we ascertained the number of late-stage ani- 
mals slaughtered for consumption, i.e. ani- 
mals slaughtered in the late stage of incu- 
bation when it is hypothesized that infected 
tissues are infectious. Between July 1980 
and June 1997, we estimated that 301 200 
(CI = [27 600-837 600]) cattle were infected 
by the BSE agent. Only 47 300 (CI = 
[3 000-135 000]) late-stage animals were 
slaughtered for consumption before the 
French specified-bovine-offal ban in June 
1996 and 1 500 (CI = [900-2 300]) between 
July 1996 and June 2000. If we restrict the 
study interval from July 1987 to June 1997, 
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Figure 2. Estimated ( — ) annual incidence of BSE infection since 1980 (A) and since 1987, using 
different scales to visualize the dynamics of the epidemic (B), assuming a reporting probability of 
5% in June 2000 and parameter p of the reporting curve = -0.6. Assuming a reporting probability 
of 10% in June 2000, the amplitude of the curve changes but its shape is unaffected. The dotted 
lines ( — ) represent 95% bootstrap confidence intervals. 



we estimated that 5 1 300 (CI = [24 300- 
84 700]) cattle were infected. In this case, 
only 6 100 (CI = [2 500-10 600]) late-stage 
animals were slaughtered for consumption 
between July 1987 and June 1996 and 1 500 
(CI = [900-2 300]) between July 1996 and 
June 2000. 

With regard to predictions, under the 
model constraint that no animals born after 
June 1997 became infected, all our fitted 
models predict that BSE cases will continue 
to be diagnosed at least until 2006 (Tab . Ill) . 



4. DISCUSSION 

The agreement of our results with those 
resulting from modelling of the British BSE 
epidemic [6, 23], with respect to the distri- 
bution of the BSE incubation period and age 
at infection, even though the BSE epidemic 
dynamics and the model used to estimate 
the parameters were different in France and 
the United Kingdom, suggests they are 
robust. As Anderson has already shown that, 
whatever the distribution used (Weibull, 
Gamma, etc.), the estimated mean was 
always in the range 4.5-5.5, and as we 
found a similar result with a Gamma distri- 
bution, we decided not to explore the 



impact of the distribution type. The esti- 
mated infection age distribution rules out 
calf milk substitutes [3] as a major source 
of BSE infection. Indeed, this distribution 
suggests that no infections occur before six 
months of age, whereas calf milk substi- 
tutes are given before this age. In addition, 
susceptibility after age one year was not 
negligible and must therefore be considered 
when designing targeted culling policies. 
Also, the age-risk function form that we 
derived resembles that postulated for vCJD 
[31], namely, no susceptibility (or expo- 
sure) in the first months (or years) of life; 
high susceptibility (or exposure) during 
young age; and then an abrupt decline in 
susceptibility (or exposure), suggesting that 
it could be a feature of both BSE and its 
human variant. 

Our results show that the proportion of 
underreporting is the most influential param- 
eter in the model and that there was substan- 
tial underreporting across the BSE epi- 
demic until rapid tests were introduced. 
Indeed, comparison of the annual incidence 
prediction with the number of observed 
cases suggests that the reporting probability 
in June 2000 was no higher than 5%. A par- 
ametric function for the time-dependent 
reporting function is a strong assumption. 
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However, in the absence of independent 
data on reporting probabilities it is not pos- 
sible to reliably estimate reporting non- 
parametrically. Nevertheless, we chose an 
extremely flexible function. Indeed, we fit- 
ted a wide range of reporting forms, ranging 
from constant reporting throughout the epi- 
demic to a reporting function with very low 
reporting rates in the early stages of the epi- 
demic and a sudden increase in the recent 
past. With reference to the reporting prob- 
ability in June 2000, the value of 5% was 
far lower than that assumed in previous 
models [19, 20] but is consistent with the 
observed rate of 20% for the period from 
July 2001 to June 2002. This value is also 
consistent with an index evaluating the 
strength of BSE clinical surveillance, namely 
the ratio of negative clinical suspicions of 
BSE reported by veterinarians to the adult 
cattle population. Between 1990 and 1999, 
the average number of negative suspicions 
was 10.7 per million adult cows per year, 
then there was a sharp increase in 2000 to 
about 80 negative suspicions per million 
adult cows per year [14]. These numbers 
can be compared to the Office International 
des Epizooties recommendations for effec- 
tive surveillance of BSE, which are to carry 
out a minimum of 100 annual investiga- 
tions of animals showing clinical signs com- 
patible with BSE per million cattle over 
30 months of age. This is ten times higher 
than the French data during the period from 
1 990 to 1 999, which indicates a lack of BSE 
surveillance in that period and thus points 
in the same direction as our results. 

The confidence intervals of the esti- 
mated total number of infections in cattle 
born between mid-1987 and mid-1997 
show less variability than estimates for the 
1980s. From 1987 onwards, the epidemic 
pattern was consistent with the control 
measures adopted to prevent the BSE epi- 
demic (Fig. 2B). Imports of British MBM 
increased strongly in the late 1980s, possi- 
bly explaining the rise in the number of 
infections between 1987 and 1990. Initial 
control measures, such as the French embargo 



on British MBM (August 1989), and the 
embargo on live cattle and the ban of 
MBM in cattle feed (September 1990), could 
explain the fall observed between 1990 and 
1992. Feed cross-contamination and recy- 
cling of infectious material within the feed 
industry may explain the rise observed 
between 1992 and 1995. It is too early to 
judge the effectiveness of measures taken in 
1994 (extension of the MBM ban to rumi- 
nant feed) and 1996 (specified-offal ban), 
based on data reported up to June 2000. 

Our estimates are different from those 
recently published by Donnelly [19], prob- 
ably owing to differences in the assump- 
tions and data used. Indeed, in her analysis, 
Donnelly considered clinical cases detected 
by passive surveillance up to June 200 1 ; she 
thus took into account the 134 BSE cases 
detected between July 2000 and June 2001 
(103 BSE cases were detected between 
1991 and June 2000). She assumed a con- 
tinuous trend in monitoring efficiency, 
despite biases induced by the introduction 
of active surveillance in mid-2000 and 
despite the fact that 179 000 cattle over 
30 months of age were destroyed between 
January and June 200 1 without being tested. 
Therefore, specific developments of the back- 
calculation method will be necessary to 
incorporate the additional cases detected 
after the outset of the targeted screening 
program. Also, these developments should 
take into account the heterogeneity of BSE 
exposure related to breeding practices. We 
could not take into account this heteroge- 
neity as almost all cases considered in this 
analysis were bred in the west of France. 
Nevertheless, the accumulation of addi- 
tional years of data, even with perfect 
reporting of all cases, would not substan- 
tially reduce the uncertainty about the early 
trends in incidence of BSE infection. In 
addition, Donnelly assumed complete noti- 
fication of BSE cases from 2001 onwards, 
whereas passive surveillance only detected 
20% of BSE cases between July 2000 and 
June 200 1 . This assumption implies a major 
underestimation of the BSE case number. 
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Comparison of estimated survival distribu- 
tions (Fig. lA) between British (-□-) and 
French (^) cattle suggests that French cat- 
tle have a shorter lifespan. Donnelly' s mod- 
els of the French BSE epidemic did not take 
this specificity into account. Indeed, she 
assumed that the "French cattle" survival 
up to 2.5 years of age was equal to that in 
the British cattle, and survival beyond this 
age was estimated from data on a single 
French cow herd (-D-). Even if survivorship 
to the age classes with most cases of clinical 
BSE (4 to 7 years) looks quite similar, the 
use of Donnelly's survival figures does not 
yield the same BSE infection numbers. 
Indeed, as almost all cattle are infected 
before one year of age and as almost all Brit- 
ish cattle survive until this age, the esti- 
mated BSE infection number will be higher 
using Donnelly's survival rates. 

With regard to predictions, all our fitted 
models predict that BSE cases will continue 
to be diagnosed at least until 2006. These 
predictions are a lower bound as they were 
obtained under the model constraint that no 
infections occurred after June 1997. How- 
ever, as of 8th January 2004, 26 BSE cases 
have already been reported in cattle born 
after that date [2]. These cases are probably 
due to cross-contamination between rumi- 
nant and not-ruminant feed, as the MBM 
ban was only extended to all animal species 
in November 2000. Hence, BSE infections 
are likely to have occurred until at least this 
date. However, the comparison of our pre- 
dictions with observed clinical cases will 
allow us to judge the effectiveness of meas- 
ures taken in 1994 (extension of the MBM 
ban to ruminant feed) and 1996 (specified- 
offal ban). Updated back-calculation will 
also be necessary to obtain updated projec- 
tions of the BSE incidence, which will 
allow the impact of potential changes to 
current control measures to be assessed. 

Between July 1980 and June 1997, we 
estimated that 301200 (CI = [27 600- 
837 600]) cattle were infected by the BSE 
agent. Only 47 300 (CI = [3 000-135 000]) 
late-stage animals were slaughtered for 



consumption before the French specified- 
bovine-offal ban in June 1996, and 1 500 
(CI = [900-2 300]) between July 1996 and 
June 2000. If we restrict the study interval 
from July 1987 to June 1997, an estimated 
51 300 (CI = [24 300-84700]) cattle were 
infected. In this case, only 6 100 (CI = 
[2 500-10 600]) late-stage animals were 
slaughtered for consumption between July 
1987 and June 1996 and 1 500 (CI = [900- 
2 300]) between July 1996 and June 2000. 

In their first study of the British BSE epi- 
demic [23], Anderson and his team esti- 
mated that 766 000 (CI = [745 000- 
799 000]) British cattle had been infected 
by BSE. Among these, only 8 000 late-stage 
animals were slaughtered before November 
1989 (date of the British specified-offal 
ban) and 43 500 in the period 1990-1995. 
In their last study [22], they estimated that 
4 000 000 animals were BSE-infected but 
the number of late-stage animals slaugh- 
tered for consumption was not specified. In 
our study, we estimated that 30 1 200 (CI = 
[27 600-837 600]) French cattle were infected 
by the BSE agent. Among these, 47 300 
(CI= [3 000-135 000]) late-stage animals 
were slaughtered for consumption before 
the French specified-offal ban in June 1996 
and 1 500 (CI = [900-2 300]) between July 
1996 and June 2000. If we restrict the study 
interval from July 1987 to June 1997, an 
estimated 51300 (CI = [24 300-84700]) 
French cattle were infected. In this case, 
only 6 100 (CI = [2 500-10 600]) late-stage 
animals were slaughtered for consumption 
between July 1987 and June 1996 and 1 500 
(CI = [900-2 300]) between July 1996 and 
June 2000. 

On the sole basis of raw estimates of late- 
stage animals in France and the United 
Kingdom, exposure of the French popula- 
tion to the BSE agent via French bovine 
products was not negligible compared to 
exposure of the French population via the 
imports of British bovine products. In addi- 
tion, the specified-offal ban on all cattle 
entering the human food chain was only 
instigated in June 1 996 in France, compared 
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to November 1989 in the United Kingdom. 
Consequently, French population exposure 
to the BSE epidemic in French cattle lasted 
until June 1996, and there was a shift of the 
temporal pattern of exposure to infectious 
animals in France and the United Kingdom. 
Nevertheless, possible differences between 
these two countries in exposure to high-risk 
bovine tissues (because of eating habits or 
risk reduction measures) should be care- 
fully considered before reaching conclu- 
sions. Further work will compare global 
exposure to the BSE agent in the French 
population with that in the British popula- 
tion, taking these differences into account. 
Finally, modelling indicates substantial 
under-reporting until active surveillance 
was introduced. This could explain why 
more cases of BSE were reported among 
cattle born after the MBM ban than before 
it. Because of this under-reporting, the 
French BSE epidemic in the late 1980s was 
completely undetected, and only the second 
wave, after 1990, was observed. Also, expo- 
sure of the French population to the BSE 
agent via French bovine products was not 
negligible compared to exposure of the 
French population via the imports of British 
bovine products. 
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